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ABSTRACT 

We present a Lagrangian-based metric theory of gravity with three adjust- 
able constants and two tensor fields, one of which is a nondynamical "flat- 
space metric" t). With a suitable cosmological model and a particular choice 
of the constants, the "Post-Newtonian limit" of the theory agrees, in the 
current epoch, with that of General Relativity (GRT); consequently our theory 
is consistent with current gravitation experiments. Because of the role of 
^ 1 , the gravitational "constant" G is time dependent and gravitational waves 
travel null geodesics of t| rather than the physical metric g. Gravitational 
waves possess six degrees of freedom. The general exact static spherically 
symmetric solution is a four parameter family and one of these solutions is 
investigated in detail. Future experimental tests of the theory are dis- 
cussed. 
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I. INTRODUCTION AND SUMMARY 

Within the past few years an elegant theoretical formalism, the "Para- 
metrized Post-Newtonian" (PPN) framework, has been developed 1 to analyze 
metric theories of gravity. The PPN framework is structured around the 
"weak gravitational fields" and low velocities of the gravitational matter 
which characterize typical solar-system tests of gravity. It classifies each 
gravitation theory as to its form "in the Post- Newtonian (PN) limit." At 
first it was hoped, and indeed seemed to be true, that the PN limit of each 
theory of gravity is unique — thus by solar-system experiments alone, one 
could, in principle, determine the "correct PN limit," which would then 
correspond to one and only one "correct theory of gravity." In addition, 
it was hoped and is hoped, that the "correct PN limit" is that of General 
Relativity (GRT) (although we try not to let this fact prejudice our investi- 
gations). To play devil's advocate, a program was initiated to attempt to 
formulate theories of gravity with the same PN limit (and hence PPN para- 
meters 1 ) as GRT. The aims of such a program are two- fold, as one can ask 
the following questions: (i) If such theories exist, how complex and con- 
trived are their formulations? (ii) Do such theories have anything in com- 
mon and in what respect do they differ from GRT outside of the PN limit? 

The first question is primarily only of aesthetic interest. But the second 
has the possibility of identifying powerful new theoretical and experimental 
tools for testing relativistic gravity — indeed that has been the case (see 
Sec. VI and Refs. 3 and 4 ). 

In this paper we present and analyze a new theory of gravity — one 
which has the same PN limit (for the current epoch) as GRT, given a suitable 
cosmological model and a particular choice of the adjustable constants. 
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Analysis of our new theory provides partial answers to questions (i) and 
(ii) above. 

A further motivation for study of this particular theory is to analyze 
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in detail the role of prior geometry in gravitation theories, a role which 
will be investigated in more general terms in another paper. 5 

To date the authors are aware of three other new metric theories which 
are candidates for sharing the property of having the same PN limit as GRT 
(candidates in the sense of contingency upon the existence of special but 
acceptable cosmological solutions and certain choices of the available adjust- 

g 

able constants). These theories are the Hellings-Nordtvedt theory, Ni's 

7 8 

theory, and the Will-Nordtvedt theory. Of these three, Ni's theory con- 
tains prior geometric elements like our own. 

A. The Lagrangian Formulation 

The equations of the theory are obtained, in the usual way, by varying 
the dynamical variables in the Lagrangian: 


L = J* d ^ x > 

(la) 

g = §(n>h) > 

(lb) 

Riem(q) = 0 , 

(lc) 


where T),h, g are second-rank symmetric tensor fields: r) is an absolute 
variable 2 (not varied in L), h is dynamical, and g is constructed algebra- 
ically from q and h. The Riemann tensor constructed out of q is denoted by 
Riem(q), and consequently Eq. (lc) states that q is a "flat-space metric." 
It is Eq. (lc), the "field equation" for q, that introduces geometrical 
structure into the theory which is independent of the matter distribution 
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— thus the "prior geometry." The gravitational Lagrangian density is denoted 

2 

by while the nongravitational Lagrangian density, * s tbe same as the 

corresponding quantity in other metric theories, (e.g., GRT), with q repre- 

A 

senting the matter fields. The "physical metric," governing the response of 
matter to gravity, is denoted by g. 

Explicitly, and g are defined by the following: 

^ = - (16,)- 1 + "vJWk- i) 1/2 ,< 2) 

V ’ (1 ' Kh)2 V A rv - (3a) 

A^ (5 “ - \ h a ) = 5 a • (3b) 

v v p 2 \x v 


Conventions and definitions for the above are the following: 

(i) Greek indices run 0-3, Latin 1-3. 

(ii) units chosen such that G = c = 1 (gravitational constant today 
and speed of light) (see Sec. Vi). 

(iii) slashes "I" and semicolons ";" denote covarient differentiation 
with respect to the flat space-metric q^g and the curved- space 
metric g^g respectively. Comma "," denotes a partial coordinate 
derivative. 


(iv) q is the determinant of q^-. 

OL 

(v) 5 is the Kronecker delta. 

' ' v 

(vi) A^ V is defined by Eq. (3b). 

(vii) indices on and h„„ only are raised and lowered with q , 

' Qp Qp 'JJ.V 

,a ,op _ , , op a ... , . 

i.e., h a = h q^g = h, and q q^^ = 5 : indices on all other 

tensors will be raised and lowered with g 

Qp 

(viii) signatures of q and g are + 2. 
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(ix) a, f , K are adjustable constants. 

Motivation for the rather ungainly expression for the metric [Eqs. (3)] 

q _ 10 

comes from an analysis of the Belinf ante-Swihart theory of gravity a 

theory which can be reformulated, at lowest order, into a metric theory 
with "effective metric" of the form of Eqs. (3). From that suggested alge- 
braic form for the metric we have constructed the present full metric theory. 

B. Summary 

Section II includes a discussion of the field equations and a calcula- 
tion of the PN limit of the theory. It is shown that there are mathemati- 
cally ten degrees of freedom in the initial value problem for h^ v (compared 

with two for g in GRT). In the PN limit there are, in general, "preferred 
|J. V 

frame effects"; 1 such effects are, however, functions of only the cosmolo- 
gical boundary values of h . By a certain choice of the cosmological model 

M' v 

one can make these effects vanish for the current epoch. We suspect that 

such time-dependent preferred- frame effects are a common property of prior 

geometric gravitation theories. At any rate, the observed absence of 

preferred- frame effects can only place upper limits on the cosmological 

boundary values of h • 

ji V 

Section III derives and discusses the equations of stellar structure 
for static, spherically symmetric stars. The equations are much more com- 
plicated than the corresponding ones in GRT (see Table i) and there is prob- 
ably no analytic solution even for a star of constant density. In addition, 
a stellar model is not uniquely specified by giving its equation of state 
and central pressure, as is the case in most other theories. The exact 
exterior, static spherically symmetric solution is obtained and is found to 


be a k- parameter family. 



Section IV includes an analysis of a special exterior spherically 
symmetric solution. For this special solution, the effective potentials 
for particles and photons are similar to the corresponding quantities in 
the Schwarzschild geometry of GRT, outside of a couple of Schwarzschild 
radii." However, the physical manifold extends only to p =. 1.5 m, 11 which 
is a "point at infinity" (not reachable in finite affine parameter by any 
geodesic) . 

There are no singularities or horizons (i.e., no black hole) in the 
physical manifold in this exact solution, but a peculiar geometrical effect 
in which the proper surface areas of concentric spheres centered on p = 0 
pass through a minimum and then increase as one moves radially inward 
(decreasing p and increasing proper time for radially falling observer). 

p 

The minimum of areas is approximately 97jtm~ and occurs near p = 2. 7 m. The 
areas then increase to infinity at p = 1.5 m, although space is not flat 
there. 

It is also found that one cannot embed the entire constant time, 
equatorial geometry in a Euclidean 3-space, but that a pseudo-Euclidean 
space is necessary for 1.5 m < p < 2. 1 m. 

Section V discusses time-dependent solutions, conservation laws, and 

12 

gravitational waves. Birkhoff's theorem does not hold in this theory, 
i.e., the exterior geometry of a spherically symmetric and asymptotically 
flat spacetime need not be static - collapsing stars can radiate monopole 
gravitational waves. The general plane gravitational wave has six physical 

3 

degrees of freedom, the maximum number possible in a metric theory of gravity. ’ 

As the theory is Lagrangian-based, conservation laws follow and one can 
construct a gravitational stress-energy complex. Appropriately defined, the 
stress energy- density of this object is positive definite for all possible 
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polarizations of plane waves. In addition there is a purely gravitational 
quantity conserved all by itself, probably of only mathematical interest. 

Section VI discusses the time dependence of the gravitational "constant" 
and further possible experimental tests of the theory. In particular, a 
search for time delays between reception of gravitational and electromagnetic 
bursts and a search for "non-GRT" type polarizations of gravitational waves 
promise to be important future experimental tests of the theory. Such tests 
would also be crucial in the theories of Refs. 6, 7, 8; and their identifi- 
cation represents an important success in our program of "devil's advocate." 


II. FIELD EQUATIONS AND POST- NEW IONIAN LIMIT 


Variation of Eq. (l) with respect to the dynamical field variable h 

V 

yields the following gravitational field equations: 

(- q) 1//2 (aO h V * a + f □ h) = - 4 jtT Q ^(- g^^dg^/dh^) , (4a) 


where 


□ tT * h“ v 


alp 


T^ = 2(- g)‘ l/2 (6 ^j/Sg^) , 


(4b) 

(4c) 


and 5 is the variational derivative. 

From the matter equations, obtained by variation of q in Eq. (l), one 

A 

can show in the usual manner (see, e.g., Ref. 13) 

1*3 = 0 . (5) 


Equation (5) is the typical "matter response equation" in metric theories. 

Contraction of Eq. (4a) with t) yields an equation for h alone, which 

M-V 

can be substituted back into Eq. (4a) to yield 
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□ h^ V = - ( 4 ir/a)(- g) 1//2 (- t])' 1 / 2 !^©^ - f(a + U)" 1 ©J q 7T t/ AV ] , (6a) 
where 

9 ^ - • (6b) 

The linearized limit of Eq. (6a) is 

□ h>* v = - (lw/«) 53 - ^ / V (f + 2 Ka) (a + Uf)‘ l ] . ( 7 ) 

Unlike metric theories without prior geometry, the four Eqs. ( 5 ) do 
not follow from the gravitational field equations; they are additional 
equations. 5 However, there is no problem of overdetermination because all 
of the 10 components of h^ V are now dynamical variables; i.e., if all of the 
essential coordinate freedom is used up in choosing a frame in which 
has a particular set of components, [usually diag(-l, 1, 1, l) ], then there is 
no coordinate freedom left to adjust the components of h . 

JJ.V 

For example, for a perfect fluid T is described by four matter 
variables once an equation of state is given (3 components of four velocity 
and energy density, for example). Thus Eqs. ( 5 ) and (6a) comprise a system 
of fourteen independent equations for the fourteen unknowns. 

We also note that all of the ten Eqs. (6a) involve second time deriva- 
tives of h . Thus in the Cauchy problem all of the h are to be regarded 

|i.V 

as dynamical variables and there are ten degrees of freedom. Once g^ has 
been constructed from q^ and h^, however, coordinate transformations can 
be performed and so there can only be six "physical" degrees of freedom. 

This is to be contrasted with GRT in which not only can four of the g^ be 
chosen arbitrarily by coordinate conditions, but also four of the field 
equations involve only first time derivatives. Thus in the corresponding 
Cauchy problem, the Einstein gravitational field has only two physical 
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degrees of freedom. 

The PPN framework of Nordtvedt, Will, and others can be used to analyze 
the predictions of all metric theories with respect to solar- system experi- 
ments (e.g., light bending, perihelion shift, gravimeter data, earth-moon 
separation, etc.). The reader is referred to Ref. 1 for a complete summary 

of the PPN framework. Briefly, this formalism involves expanding the metric, 

II4. 

in the manner of Chandrasekhar, in the small dimensionless quantities which 
occur in the solar system stress energy tensor, e. g. , 

v 2 ~ U ~ (p/p) ~ n ~ 0(e 2 ) M lo -7 , (8) 

p 

where v is the squared velocity of a typical fluid element, U is the 
Newtonian potential, P/p is the pressure divided by energy density (specific 
pressure) and n is the specific internal energy. It is found that, in a 
particular coordinate gauge, and for most metric theories — including ours 
— there are only nine different functionals which can occur in the metric 
at PN order and only nine independent parameters multiplying these functionals. 
Almost all twentieth century gravitation experiments to date can be summa- 
rized by their constraints on these nine parameters, the "PPN parameters." 

We now calculate in our theory the PN limit, which will involve a 
perturbation solution of Eq. (6a). For calculational ease we assume a 

coordinate system in which q^ takes on Minkowski values. Before we begin, 

15 

a crucial point must be recognized. The metric has the form 

8 q£ = V + ' 

and we know that far away from the solar system there is some coordinate 
system in which g^ takes on Minkowski values. However, this coordinate 
system will, in general, not be the same frame in which q^ takes on 


8 



Minkowski values; there is no a priori reason why the boundary values of 

h should be zero in this coordinate system. Thus in solving Eq. (6a) 

|J.V 

we are not at liberty to set equal to zero for all time the "arbitrary 

constant" which may be added to h ; this complicates considerably the con- 

p-v 

struction of the PN limit of our theory. However, we feel that this com- 
plication and its origin are of sufficient educational value to warrant a 


detailed discussion. 

Denote the nearly constant boundary values of h^ v by (u>^ v can 

only change on a cosmological time scale by definition) and the part tied 

* 

directly to the solar system by h : i.e., 

r " 



+ U 


HV 


(9) 


Now use the six-parameter invariance group of the Minkowski metric to 
pick a coordinate system in which u is diagonal, reducing w to four 
components. Without justification, but for simplicity, we now assume that 
the three spatial components of w are equal. Such an assumption does not 

pi V 

effect the qualitative conclusions of this section. Further assume that 


!%,! « 1 > ( 10 > 

although does not have to be as small as the 0(e) indicated in Eq. (8). 
Equation (10) will turn out to be an assumption consistent with the ultimate 
experimental limits on the w . 

r* " 

Next expand Eqs. (3a) and (3b) in a power series in h : 

o = ti — 2Khn + h + K^h ti — 2Khh +T-h hv+... . (if) 

®|iV >V Vv l^V '(iV |IV 4 (IT 

When Eq. (9) is substituted into Eq. (ll) one obtains 
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(12a) 


g 00 °0 + E o h oo 


* p p * * 

- F 0 h - K h - 2Kh 00 h - ^ h Q0 


g ij = DS ij + Eh ij + F5 ij h “ 2Kh 


’Ok ^Ok * 


♦ 1 »;/ - 


(12b) 

(12c) 


where all of the constants appearing in Eqs. (12) have the form: 


Dq = 1 + 0(co), etc., and are given explicitly to 0 (cj ) in Appendix A, along 
with other constants appearing below. Using Eqs. (12) and a perfect fluid 
for the matter stress-energy tensor, one obtains from Eq. (6a) 


P 


□h*^ v = _ (U.ji/a) l“ 1 pv°v^(l + + Igh + I 3 v 2 )|(l - 2Kto) 5^^ 

+ l^ v + 1 s V P + M ^ Vw op + 

+ 1 & V* v g - “V/, + W'Vj • 

In Eq. (13) I, 1^ Ig, Ig, M, N are all functions of a, f, K, (see 

Appendix A) and 


(13) 


w = 3a) n - w Q0 , 


v a = dx a /dt , 


(14a) 

(14b) 


P = proper mass-energy density measured 
in the rest- frame of the fluid. 


(14c) 


To simplify an already complex presentation, we have omitted the pressure 
from the perfect fluid stress energy tensor and included the internal energy 
in the total proper energy density p. (Such terms are not omitted in quot- 
ing the final PPN parameters. ) We now write 

h *pv = (l) h *pv + (2) h *nv + _ ? (15) 
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in a perturbation expansion and obtain (see Appendix A for notation) 


y 2 (l) h *00 = _ _ 2 Kw) + L - u> 0 (| + M)j = - hnp C Q 

, (16a) 

V 2 ^h* 1 ^ = - ^rtpx(Mco 0 - L) s lj £ - CjB 1 ^ , 

(16b) 

V 2 ^h*° k = - lj-riTp|v k (l - 2Kw) + co^v^ J = - 4itp c 2 vk > 

(16c) 

v 2 ( 2 V 00 - - ws 0 (1) h*°° + Sl (1 V + B 0 v 2 ) + CDh^ 

, (l6d) 

y 2 (2) h *ij = _ + 8 lj (R 1 ^ 1 ^h* 00 + R^h* + 

2)1 + (l) h *U 
; | ,00 

where 


r =, (al) _1 . 

(17) 

Solutions of the equations are 


^V 00 - 0 0 M , 

(18a) 

(Dh* 1J = 5‘V , 

(18b) 

(1) h*° k - c 2 v k , 

(18c) 

(2)h * 00 - T l S 0 C 0 + S 1 (3C 1 - VK + TB 0»1 + c 0 x ,00 ’ 

= tR 0 4 J + T S iJ | R 1 C 0 + R 2 (3C 1 - c 0 >) *2 

(I8d) 

+ TBjS 11 ®! + C lS Ij X ;00 , 

( 18e) 

where we have defined the five "potentials" U, V^, (J’p $ 2 , > 

and the 

"superpotential" x aa follows: 


U(x, t) = J* p(x', t) |x - x' | _1 d 3 x' , 

(19a) 

V fc (x,t) = J p(x', t)|x - x' J " 1 v k d 3 x' , 

(19b) 

$ x (x, t) = J* p(x’, t) v 2 |x - x’ I " 1 d 3 x' , 

(19c) 


, ( 16e) 
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(I9d) 


® 2 (x, t) 5 J p(x', t) |x - x' I" 1 U(x', t) d 3 x' 
<J lj (x, t) s J* p(x',t)|x - x' I" 1 vV d 3 x' 
v 2 x = U . 


(I9e) 
( 19f ) 


Using Eqs. (12) and our solutions, Eqs. (18), we now compute the metric: 


8 oo - - »o + h v + k 2 ° 2 + Vz + Vl + *1*00 ’ 

S ij = 6 ij (D + K 5 U) ' 

®0k = ” HC 2 V k * 


(20a) 

(20b) 

(20c) 


Notice that the metric does not approach the standard Minkowski tensor far 
away from the solar system (when the potentials U, <j>p $> 2 ' X ^ 0) because 
of the leading constants D Q and We must therefore make a "scaling" 

transformation: 


t = Do" 1 / 2 t , (21a) 

x = x . (21b) 


In the tensor transformation law for the metric 


V® * 8 °e® SJ* 


dx a Bx^ 


8q 8 [U( ~’ t}> •l ( 5' t) ' 

dx 


] 


dx a 

S? 1 dx v 


( 22 ) 


we also need to express the potentials as functions of the new (barred) 
coordinates. An example of the procedure is the following: since p is a 

scalar 


p(x, t) = p(x, t) , ( 23a ) 

U(x, t) = J* p(x',t)|x - x' I ' 1 d 3 x' = J p(x',t)|x - x' I " 1 d 3 x' , 

= D _1 J p(x', t) |x - x' f 1 d 3 ^ = D' 1 U(x, t) . (23b) 
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In a similar manner one finds 


<D 2 (x, t) = 

_ 2 

D $ 2 (x,t) , 

(23c) 

$ 1 (x, t) = 

DqD" 2 ^^,!) , 

(23d) 

< 

7? 

rt 

II 

D 0 1 / 2 D" 3 / 2 V k (x, t) , 

(23e) 

X 

o 

o 

II 

2 V - - • 

w r\r\ 

(23f) 


Making the transformation indicated in Eqs. (22) and (23) and then dropping 


the bars, g becomes 
- 1-1 


- 1-2 2 
c —tT TT 4 — 


-1-2. 


g 00 = - 1 + V D K 1 U + D 0 ° K 2 u + D 0 ° V2 + D Vl + U K l\oO ' 


+ D"TL 


*lj - » t J (1 + D ' 2 k 5 U) ' 

*0k ’ - HC 2 D ‘ 2v k • 

A final coordinate transformation must be made to remove the 0Q term from 

e and reduce the metric to "standard PPN form." However, additional 

®00 

transformations of the form of Eqs. (23) are now negligible corrections and 

no distinction need be made between functions of new and old coordinates. 

-2 

The result of the final transformation, t -> t +l/2 D 

*00’ *00- V^OO ’ (25a) 

8 ij " g ij ' 

*0k ’ *ok + ¥ K i D ‘ 2(v k - V > 


(25b) 

(25c) 


where W, is a new potential defined by 


W, = J* p[v • (x - x')] |x - x' I" 1 (x - x' ) k d 0 x' . 


rl ^ • 


(26) 


(24a) 

(24b) 

(24c) 
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We now demand the proper Newtonian limit, i.e., 


which requires 

l^D - 1 D -1 = 2 today ; (27) 

(a consequence of our choosing units in which the gravitational constant is 
unity today). Equation (27) expresses a (constraint between the three adjust- 
able constants a, f, and K for a given set of u> . Comparing Eqs. (2l0-(25) 
with the definitions of the PPN parameters and using Eq. (27) to simplify, 
one finds 


y = ■ \ D' 2 !^ = 7 (a, f, K) + 0(w) , 

(28a) 

P = - \ D 0 -V 2k 2 = P(a,f,K) + 0(u) , 

(28b) 

5i - c 2 “ S s = ^ = a 3 = 0 ' 

(28c) 

a x = 2HC 2 D' 2 - by - k = 0(u) , 

( 28d) 

a 2 = DqD -1 - 1 = 0 (u) . 

(28e) 

where y and p are defined implicitly by the relations 


a = (27 + 2)" 1 , 

(29a) 

f = (lop + 67 P - iy 2 - Qy - 6)[2 (y + 1)(37 - 

5 - bfi) 2 ]' 1 .(29b) 


In GRT, 7 = p = 1 and the other seven parameters vanish. In our theory it 
is clear that the two adjustable constants, a and f, may be so chosen to 
give any value to y and p. For example, if the u are all zero, one can 
satisfy Eq. (27) and have y = P = 1 with the choice 

(a, f, K) = (p - -j^) • ( 30) 


lb 



leads to non- 
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It has been shown that the nonvanishing of Q^, or 

invariance of the functional form of the metric of Eqs. ( 2k)-(25 ) under 
post-Galilean transformations (curved-space versions of Lorentz trans- 
formations). New terms, involving the velocity of the Lorentz boost with 
respect to the current "preferred frame" and multiplied by combinations of 
a v Qtg, a y appear in the metric. Nordtvedt and Will have calculated the 
experimental consequences of the resulting "preferred- frame effects" and 
find that they lead to periodic anomalies in such phenomena as the solid 
earth tides, secular perihelion shifts, etc. The reader is referred to 
their paper for further details and we quote here only the current experi- 
mental limits on 0!^ and OC^: 


0 !^ — 0. 1 , 


a 2 ^ 0.02 . 


(31a) 

(31b) 


We have calculated explicitly the quite complicated functions (^(w^), 

a (to ) and have examined their numerical values over a large range of con- 
2 V |iv 

stants a and f (consistent with the experimental limits on j and £). We 
find that the experimental constraints indicated in Eqs. (31) require appro- 
ximately 


|co Q | + 1^1 £ .015 . (32) 

Even if we had not made the simplifying assumptions about the form of oj^, 
its individual elements presumably would still be required to satisfy roughly 
the constraint of Eq. (32). 

Since the w are cosmological boundary values of h , one must solve 
(iV 

the cosmological problem for a particular cosmological model to obtain the 
theoretical values of the w . Because of the absolute nature of it 
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should be possible to construct cosmologies such that, during the current 
epoch, the curved and flat-space metrics approach Minkowski form, far from 
the solar system, in the same coordinate system. Such a cosmology would 
guarantee that the u> vanish at present, although a time dependent cos- 
mology would certainly cause nonzero values of to occur over cosmolo- 
gical time scales. Indeed, perliminary results from a cosmological solution 
possible to make all of the u^ v arbitrarily small for the current 
indicate that it is/epoch — and still have a reasonable cosmological model. 

Thus, a consistent solution exists for which the PN limit of our theory is 
arbitrarily close to that of GRT in the current epoch. 

Further details regarding the time dependence of the are given in 
Sec. VI. 


III. THE GENERAL STATIC SPHERICALLY SYMMETRIC 
SOLUTION AND EQUATIONS OF STELLAR STRUCTURE 
A. The General Exterior Static Spherically Symmetric Solution 


Before writing down the equations of stellar structure for a static 
spherically symmetric star, let us construct the general static spherically 
symmetric exterior solution (which must then be joined onto the solution 
inside the star). 

First of all, choose a coordinate system in which 


’l 


pv 



1 


2 

r 



( 33 ) 


The most general form of 
symmetry requirements is 



in this coordinate system which satisfies the 
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0 


0 


A(r) 


r^sin^O / ( r ) 


The homogeneous field equations for h are simply 

Qg, (IV „ 

n h = 0 . 


The solutions to Eqs. (35) which are well behaved at infinity are £ 


/ V r , 

( Y 


- 2 V r , 

Jx - 2a 3 /r" 


r 2 (a 0 /r + a ,/r 3 ) 


r 2 sin 2 G(a c /r + a ,/r 3 ) ,(36) 


where a , a 0 , a_, and a, are arbitrary constants. We remind the reader that 
the r coordinate in Eq. (36) has, at this point, no interpretation other than 
its relation to the group — theoretically defined assumption of spherical 
symmetry. Construction of g from h is purely algebraic [see Eqs. (3)], 

t-L V |J-V 

and the details will not be given here. Since h^ v has off-diagonal terms, 
so will g . However, having obtained g^ v , we can make the coordinate trans- 


formation 


t -*■ t + 


- . 


which then diagonalizes the metric, and we finally obtain 


„ _,2 2 [ a l ( /. 1 a 2 a 3\ 2 

(1 - Kh) 7 — T ' ( 1 ' 2 T Ti ’ 

r r ' J 


(38a) 
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Equations (38) for the metric indicate a 4-parameter family of solutions for 
the general static spherically symmetric exterior metric. One can convince 
himself that all four of the parameters are physical (not removable by coordi- 
nate transformations) by transforming to curvature coordinates and verifying 

22 

that four arbitrary parameters remain. In Sec. IV we will investigate 
more closely a particular member of the 4-parameter family. 

B. Stellar Models 

We idealize a star as a spherically symmetric, static mass of perfect 
fluid and assume a temperature- independent equation of state 

P = p(p) , (^°) 

where p is the pressure and p the energy density. We work in the coordinate 

23 

system in which tj has the form of Eq. (33). For mathematical simplicity 
M* v 

we seek solutions for h which are diagonal, i.e., with p(r) = 0 in Eq. (34). 
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Such solutions represent a subclass of all possible solutions and corres- 

2b 

ponds to the condition a^ = 0 in the exterior metric [cf. Eqs. (36) and (38)]. 
The metric now has the form 



Equations (5) and (6a) together with Eq. (40) are the necessary set for 
computing the structure of our stellar model. With the usual fluid stress 
energy tensor 

T°^ = (p + p) u a u^ + pg°^ , (^2) 


one finds that the only nonvacuous equation resulting from Eqs. (5) is 

dp/dr = -|(p + p) [ 2K( 1 - Kh)" 1 dh/dr + (1 + -gCp)" 1 dcp/dr] . (^3) 

Using the Christoffel symbols for q, one finds that Eqs. (6a) yield 
the following: 

! 

- 3p) + (1 - Kh) |p( 1 - \ t)” 1 + 2p(l - A) -1 

(44a) 

V 2 cp = (f/a^ 2 *! - (lnt/ a ) r [ 2K ( 3 P - P) + (! + -§«p) 1 C 1 “ Kh )p) » 

V 2 ^ = 4(\|r - A) r' 2 - (f/a) - (4*/a)r[2K(p - 3p) + (1 - \ t) _1 (l - Kh)p] , (44c) 


v 2 ** = - 4itr(a + 4f) " 8 

- p(l + -^p) 1 ] , 


where 
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V 2 5 d 2 /dr 2 + 2r -1 d/dr , ( 45a ) 

ra(l- Kh) 3 (1 + ^) _1 (1 - MO' 1 (1 - 4 A)' 2 • (45b) 

Equation (44a) follows from taking the trace (with respect to rj) of Eq. (6a). 
Equations (44b) and (44c) are the 0-0 and r-r components respectively of 
Eq. (6a). Altogether, Eqs. (40), (43)- (44) are five highly nonlinear coupled 
equations for the five unknowns p, p, Cp, A, and i|/. . Linear combinations of 
Eqs. (44) can be taken to yield 

V 2 (* - *) .[-|(r/a)(l - Kh)(l -- I *)' 1 (1 --!*)' 1 + 6r- 2 ](* - A) , (1*6) 

which is an equation we will later discuss. 

Outside of the star the physically acceptable solutions to the homo- 
geneous forms of Eqs. (44) are [cf. Eq. (36)] 


cp = a x /r , 

(47a) 

i = a 2 / r ' ’ 

(47b) 

\ = a 2 /r + a 3 /r 3 . 

(47c) 


The constants a p a 2 , and a 3 are to be determined by matching conditions 
at the surface of the star. The general procedure in constructing stellar 
models is to choose various central values for the variables, integrate the 
equations outward from the center until the pressure vanishes, and thus 
establish the surface of the star. Various boundary conditions must typi- 
cally be satisfied, but in the case of GRT, for example, the conditions can 
be satisfied in a trivial manner without multiple trial integrations. The 
situation here, as we shall see, is vastly more complicated. 

As long as the denominators do not vanish (see discussion below), 
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Eqs. (44) are regular at the stellar surface and hence require that rep, r if, 
rA and their first derivatives be continuous across the surface. Using Eqs. 
(47) and denoting quantities evaluated at the surface by a subscript s, one 
obtains the six matching conditions: 

cp g = a^R, (cp ?r ) s = - 3l /R 2 , (48a) 

^s = a 2^ R ' ’ ^,r^s = " 3 2/ r2 + 6a 3/ R ^ * ^ 8b ) 

A g = a 2 /R + a 3 /R 3 , (A^) g = - a 2 /R 2 - 3a 3 /R U , (48c) 

where r = R is the surface of the star. 

What are the appropriate central quantities to be specified? Suppose 
we regard (\| r - A), and cp as the three independent gravitational potentials. 
Then a possible but nonunique solution to Eq. (46) is - A = 0 everywhere, 
corresponding to considering r an isotropic radial coordinate. However, 
forgetting this special case for the moment, the regular solution of Eq. 

(46) near the origin is 

2 

\|r - A ~ const, r 

Thus one central condition to be specified is 

[(*-*)/r 2 ] c , 

where we denote by c quantities at the center, analogously to the quantities 
at the surface discussed above. The equations for h and cp are regular at 
the origin as long as the potentials are sufficiently small and therefore, 
in analogy with the corresponding electrostatic equations, the derivatives, 
at the center, of cp and A must vanish. However, the central values of the 
potentials themselves must be specified, and hence the two other central 
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parameters are <p c and A c - Thus in general we have six parameters to adjust, 

O 

e.g., a , a , a , cp , A , [(\|r - A)/r ] in order to satisfy the six matching 

constraints given in Eqs. (48), for a given equation of state and central 

pressure. One way of viewing the boundary conditions is that ep c , A £ , 

[ _ A)/r 2 ] c must be so chosen as to match onto a regular exterior solution 

at the star' s surface — such a two-point boundary value problem in general 

has a discrete set of solutions, i.e., for a given p c and equation of state 

there may be no t cp ^ A, [ (i)r - A)/r ] I such that there is a solution, or 
lee c 

there may be many different sets. Thus the central pressure and equation 
of state do not uniquely specify the stellar model in general. However, we 
do know that for a weakly gravitating star (p c /<P> « 1, cp* A, « l). 
Equations (MO become linear and do indeed have unique and well behaved solu- 
tions for each central pressure (Newtonian, and post-Newtonian regimes, see 
Sec. II). However, we can expect that as the models become more and more 
relativistic, a point is reached where each p^ and equation of state branches 
into a discrete spectrum of stellar models. 

If one trys as a solution to Eq. (46) ^ = A, then a more convenient 
form of the boundary condition is 

[ cp/(rcp )] = - 1 , 0+9 a ) 

> r S 

[(3A + rA^ r )/3A] g = f . (49b) 

One then adjust A^ and cp^ to satisfy Eqs. (49) and defines a^ and a 2 
(a 3 = 0) by 

a n = Rep , (50a) 

J. s 

a 2 = R(3A + RA >r ) s • (50b) 
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If \|r / A, then the proper constraints are 


[cp/rc^ r ] g = - 1 , (51a) 

[ (3A + r^ r )/(2X + i|r)] g = -| , (51b) 

R[ (i r - A )/ (A - Mf)] = 3 , (51c) 

> r 1 r S 

2 

and one adjusts A , cp and [ (\)r - A)/r ] to satisfy these three constraints; 
c c c 

defining and a^ as in Eqs. (50), and 

a 3 * I R X - * S ) ' (52) 

As far as the exterior metric is concerned, all of the information 
about the stellar model is contained in the parameters a^, a 2 , and a^ (and 
in the general case). Each different set of values for these constants 
corresponds to a different mass and radius of the star. Indeed, the total 
mass-energy of the star ("gravitating mass") as determined by g Q0 and using 
Eqs. (Ul) and (^7) is 

m ~ ^ a l + K ( 3a 2 “ a i) > ( 53 ) 

( a ^ and a^ determined by matching conditions at the surface). It is dif- 
ficult to say what each parameter corresponds to physically (in terms of 
integrals over the source, etc. ) because of the complexity of the inhomo- 
geneous equations [cf. Eqs. (Ml-)]. The only definite statement is that the 
particular combination of a^^ and a given in Eq. (53) corresponds to the 
total mass. 

A further interesting fact is that, for a given choice of a, f, K, the 
PPN parameters 7 and p — as determined by a l/r expansion of the isotropic 
version of the metric — are functions of a^^ and and in general are not 
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equal to their values as determined in the PN limit. (This situation is 

25 

also true in the Dicke-Brans- Jordan theory. ) Only in the case of a weakly 
gravitating star can one be sure that the two different determinations of y 
and p will agree approximately (to within PN precision). In GRT, on the 
other hand, expansion of the Schwarzschild metric gives y = P = 1 regardless 
of stellar model, and in agreement with the y and P as determined in the PN 
limit of the theory. 

Table I gives a comparison between our stellar- structure equations 
and those of GRT. 


IV. ANALYSIS OF AN EXACT EXTERIOR SOLUTION 

i 

A. Hie Metric 

As pointed out in the last section, the general exterior metric of a 
static spherically symmetric spacetime is a 4-parameter family [cf. Eq. (38)]. 
Let us analyze a member of that family. First of all, for simplification, 
we choose a^ = a^ = 0, which puts the metric of Eq. (38) in isotropic form. 
Next, using Eq. (53) as a definition of the mass m, we choose a a 2 , and K 
such that a l/r expansion of the metric indicates that the PPN parameters 
7 and P are both unity (see Sec. II). In other words, choose aj^ a 2 , and K 
such that 


8 00 

- 1 + 2m/p - 2(m/p) 2 + 0(p 3 ) , 

(54a) 

8 ij = 

- 5 i j(l + 2m/ p) + 0 (p" 2 ) > 

(54b) 

which requires 

a 1 /m = 1 , 

(55a) 


a 2 /m = 3 , 

(55b) 


K = 1/16 . 

(55c) 
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It is interesting to note that the value for K given in Eq. (55c) is the same 
value required for 7 = 8 = 1 in the weak-field PN expansion [cf. Sec. II and 
Eq. (30)]. Using Eqs. (55) and Eq. (38), one can now write the line element 



The line element given in Eq. (56) is the simplist static spherically sym- 
metric metric which yields the same light bending and perihelion shift (viz. 

7 = 8 = 1) as in GRT. (Note that the value of g Q0 is identical to the correspond- 

pC . 

ing term in the isotropic form of the GRT Schwarzschild geometry. ) 

B. Geodesic Completeness and Radial Geodesics 

A glance at Eq. (56) reveals that p = 1.5 m is an infinite proper radial 
distance away from any p > 1.5 m. To investigate whether this point is 
removed from the physical manifold we need to look at null and timelike 
geodesics. Consider equatorial orbits (no loss of generality with spheri- 
cal symmetry) and consider the first integrals of the motion for particles 
and photons: 


A, - - 1 - * + 8pp< uP > + , 

W P P > 2 + s°°(V 2 ♦ «*V = 0 ' 


(57b) 


where u a = dx a /dT for particles and p a = dx^/d A (with A the affine para- 
meter) for photons. It is well known (see, e.g., Ref. 27) that for a metric 
of the form of Eq. (26), u Q , u^ and (P are all constants of the motion, 
which we shall denote by - E, L, and i, respectively. Physically, these 
constants are energy per unit rest mass, angular momentum per unit rest mass, 
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and impact parameter respectively. 

Using the above, Eq. (57a) can be written as 


U p = dp/dT = (p - |)" 2 (p + -|)Cp - f)[E 2 - r 2 (L,p)] , 

(58a) 

where 


! 

p = q fvcij Li — L/m, etc* y 

(58b) 

and 


r(L,p) * (5 + 1)' 1 ! (p - |) 2 + (I/p) 2 0p - f) V /2 • 

(58c) 


The function r plays the role of an effective potential, which we shall 
discuss later. Equation (57b) can be written as 


where 


(dp/dt ) 2 





7 


s p(p + -|)(P - -§) 1 • 


(59a) 


(59b) 


Consider first radial geodesics (L = l = 0). Then Eq. (58a) indicates clearly 
that "p = 3/2 is an infinite proper time away from timelike geodesics. If 
one then uses the fact that P Q = g Q0 dt/dA = constant for the null radial 
geodesics together with Eq. (59a), then it is also easy to show that 
p = 3/2 is an infinite affine parameter distance away for null radial geo- 
desics. Equations (58) and (59) indicate that nonradial geodesics between 
any two values of "p take even longer proper time and affine parameter than 
do radial geodesics. Thus we have shown that p = 3/2 is really unreachable 
by particles and photons; in particular, the manifold covered by our coordi- 
nate system is maximal . 28 Since one can also show that there are no singu- 

28 

larities for "p - 3 / 2 , our manifold is geodesically complete. 

For the special case of radial geodesics, we integrate Eq. (58a) to 
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yield 


- + 


- b” 1 fn 


2b (X + 1) + 2c (j) - |) 
— - 

P - O 


- d X + d 


-3/2 sin -l/2dp^_l^ 


~2 , 


V 


(60a) 


2E 


+ const. 


for | < 


E < 1 




where 


b = (Us 2 - l) 1 / 2 , 

c s 2E 2 - 1 , 

d s 1 - E 2 , 


X *.[(1 + E 2 ) p - cCp 2 + \)} l f 2 • 


(60b) 

(60c) 

(60d) 

(60e) 


We will not be interested in analytic solutions for values of E other than 
those indicated in Eq. (60). To obtain the functional relationship between 
coordinate time t and p for l/2 < E < 1, add to Eq. (60a) a factor of b 
multiplying the log term and a factor of (3 - 2E ) multiplying the inverse 
sine term. 

For radial photon geodesics, Eq. (59a) can be integrated to yield 

3 

*t = ± Cp + 2fn|"p - §|) + const. . (61) 


Figure 1 illustrates a few of the radial geodesics for photons and particles, 
the latter released from rest at p> = 10 and "p = 5. It is interesting to note 
that the analogous metric in GRT is geodesically incomplete: p = l/2 can be 

reached in finite proper time, but requires infinite coordinate time. 

It can be shown, from analysis of the metric, that another complete 
universe exists for l/2 ^ ‘p ^ 3/2. However, if we assume the geometry to be 
produced by a star which originated in our universe, then its surface lies 
outside p = 3/2. In the following we consider only the region p > 3/2. 



B. Proper Surface Areas and Embedding Diagrams 




There are some curious geometrical effects in our manifold, not to be 
found in the Schwarzschild geometry of GRT. The proper surface area of a 
sphere described by p = const, is 


A . 2 _2 1\2 3x-2 

A = kirn p (p - g) (p - g) 


(62) 


>• A plot of this area is given in Fig. 2, in which the abcissa is marked off 
not only by p but also by the proper time as measured by a radially falling 

observer. As can be seen in the figure, the observer sees the sequence of 

2 

surface areas pass through a minimum, (h9/b + 5/6) at 

~p = 3/2 + I/2/6, and then increase without bound as p = 3/2 is approached. 

Another interesting feature arises when we examine the intrinsic geo- 
metry of the 2-surface: t = const., © = it/2 by the use of an embedding 

diagram. By equating the two-dimensional metric 


.. 2 1x2 , 3\-2 ,, 2 2,.2v 

ds = (p - g) (15 ~ g) (dp + p d© ) 

to the metric of a surface of revolution in a Euclidean 3-space 
ds^ = dz^ + dr^ + r^d<p^ =[(dz/dr)^ + 1] dr^ + r^dcp^ , 


(63a) 


(63b) 


one can visualize the geometry of Eq. (63a). If we can find z(r), or more 
easily z(p) and r(p), then the line element of Eq. (63b) can be drawn. 
Clearly 


r - p(p - -|)(p - ^r 1 


(eh) 


The function z(p) is the solution of the equation 

/dz\2 - /d s \ 2 (*lY 

\dp/ \dpj ydpj 


(p - j ) 2 [ P 2 - gp + ff 

CP - f ) 2 " 


[ P - 5p + f l 

"l (P-l 2 


(65a) 
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or 

dz/<to - o 1 / 2 (2p 2 - 5p + |) l/2 (p - |)‘ 2 • (65b) 

The right-hand side of Eq. (65b) becomes complex at 

p = -i( 5 + Z 13) 2.1 or r *» h.8 . (66) 

This indicates that for 1.5 < p < 2. 1 we will have to embed in a pseudo- 
Euclidean space, i.e., 

ds 2 = - dz 2 + dr 2 + r 2 d<p 2 . (67) 

^ The embedding diagram is given in Fig. 3 and includes both the Euclidean 

part and the pseudo-Euclidean part. The surface is obtained by rotating 
the curve about the z or iz axis. 

C. Particle and Photon Orbits 

Analysis of orbits is facilitated by use of the effective potential. 
Equations (58c) and (59b) give the effective potentials for massive particles 
and photons. For a given value of L, the particle is allowed only in those 
regions for which r(L,p) ^ £. For photons, y 2 acts as an "inverse" effec- 
^ tive potential; photons are allowed only in regions for which y - l. Figures 

and 5 illustrate the effective potentials for particles and photons, respec- 
tively, with the dots in Fig. k indicating extrema of the potential (circular 
orbits). The closest stable circular orbit for particles occurs for L ~ 3.88 
at p ~ 7. For particles with larger L, the circular orbits with p < 7 are 
unstable and those with p > 7 are stable. The circular photon orbit occurs 
at p = 1.5 + /3 or r ~ 5. This can be compared with the corresponding value 

of 7 = 3 in GRT. 
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I 

V. GRAVITATIONAL WAVES AND CONSERVATION LAWS 
A. Monopole Waves 

In the full theory (no linearized approximation) the homogeneous field 
equations are, as indicated previously, 


OP, |iv 
T) h 


a |p 


= 0 


( 68 ) 


and gravitational waves travel geodesics of q rather than g. The impli- 
cation of this last fact will be explored later. The simplicity of the 
vacuum field equations [cf. Eq. (68)] is of great help in constructing 
solutions. 

Consider a time-dependent spherically symmetric solution to Eq. (68), 
for example 


00 


= r 


■ 1 e i w(r— t) 



-1 iw(r-t) 
r e v 


(69a) 

(69b) 


The Riemann tensor constructed from the resulting time-dependent spherically 
symmetric metric is itself time dependent. From this we conclude the pre- 
sence of physical monopole waves; thus there is no analogue of Birkhoff's 
theorem 12 in this theory. The existence of such solutions in our theory 
and the accompanying monopole radiation complicate the problem of the spheri 
cal collapse of a star. As will be shown below, there are other "non-GRT" 
type gravitation- wave modes in addition to the monopole waves. 


B. Linearized Theory and Plane Gravitational Waves 

In analyzing weak gravitational waves, one should restrict one's atten- 
tion to the form and behavior of the Riemann tensor, not only because it is 
gauge invariant (under infinitesimal coordinate transformations) but also 
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because it is that feature of the gravitational wave which interacts directly 

with test bodies. Work in a coordinate system in which q is Minkowskian 

and h is small (small deviations from flat space). Then 

pv 


z = ti + h - 2KhTi 
5 pv 'n v M- v 'liv 


+ 0(h 2 ) = 


‘p.v 


+ h' + 0(h 2 ) 

[IV 


(70) 


and 


d - — fh 1 
0^75 2 V 05, P7 


+ h 3r,os 


- h 


ay, Os 


h P5,ay ) 


(71) 


Furthermore, restrict one's attention to those solutions of Eq. (68) which 
represent plane waves travelling in the z direction, i.e.. 


h' 

\iv 


= A e 
M.V 


Lk(z- t) 


(72) 


where A is a constant amplitude and k a wave number. To analyze the 

(JLV 

decomposition of R^ i nto independent "wave modes" in as invariant a manner 
as possible, one should investigate the transformation properties of R^^ 
under those Lorentz transformations which leave the wave direction fixed. 

With such transformations in mind one selects a new basis in which the com- 


ponents of R g are to be computed — the qua si- orthonormal tetrad basis 
(see, e.g., Ref. 29 for a complete discussion of the "tetrad formalism"). 


k = 2 -1 / 2 (-l,0,0,l) , 

(73a) 

l = 2 -1 / 2 (l,0,0,l) , 

(73b) 

m = 2 _1 / 2 (0,l,i,0) , 

(73c) 

m = 2 _1 / 2 (0, l,-i,0) . 

(73d) 


Note that one of the "tetrad legs" points along the direction of the wave. 
In such a basis the components of the Riemann tensor are 


^nkmf 


5 


n ‘W* 8 


etc. 


(74) 
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Using Eqs. (7l)-(74) one finds that the only nonvanishing components of the 
Riemann tensor are those with two £'s — thus there are six possible degrees 
of freedom. Since there are no restrictions on the Riemann tensor once 

Eqs. (68) are satisfied, all six tetrad components will in general be non- 

! 

vanishing and our theory thus has six independent gravitational wave modes . 

In GRT, as a contrast, the field equations V- ° imply vanishing of 

r , R , R , — , and R„ .—so that there are only two degrees of freedom 
£k£k’ £k£m’ £k£nr £m£m 

— those represented by R fm£m and its complex conjugate R f — . 

The reader is referred to Refs. 3 and 4 for details of the transformation 
properties of the objects indicated in Eq. (74). Here we quote only the 
results: We denote the six wave modes by 1^, Tj, <I> 22 and in 

terms of the tetrad components of the Riemann tensor and "electric" coordi- 
nate components of the Riemann tensor (those which are directly physically 
measurable) these are 


¥ 2 "6 R £k£k 

T 3 = " 2 R £k£m 
*3 * - \ R £k£m 

\T/ — — P = 

4 ” £m£m 

\ S - R £m£m = 

. -in __ 

$ 22 ” 2 £m£m " 


6 R tztz * 


2^ R txtz ~ ^ R tytz^ ’ 


■ + "W ’ 


_ 

R + 

2iR„ „ 

tyty 

tztz 

txty 

1 

R 

2iR„ ^ 

tyty “ 

tztz 

txty 

R txtx 

“ R tyty 

# 


(75a) 
(75b) 
(75c) 
( 75d) 

(75e) 
( 75f ) 


The presence or absence of a T 2 component in a gravitational wave is Lorentz 
invariant. If Tg is absent in a particular wave, the presence or absence 

of if (or \jr ) in that wave is also Lorentz invariant. As outlined in Refs. 

3 3 
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3 and 4, if either y 0 or y is present in a wave (in many theories they are 
always absent, but not ours), then it is impossible to decompose the wave 
into states of definite helicity (spin) in a Lorentz invariant manner: what 

one observer identifies as '’pure spin 0" another observer will identify as 
"pure spin 0" plus "pure spin 1," etc. . Only waves containing only $ 22 , 

and can be decomposed into pure spins: spin 0 and spin 2. In general, 

then, there is no unique spin decomposition of waves in our theory and it is 
of class II (see Refs. 3 and 4 for a complete discussion of the "classifica- 
tion scheme" ) . The physical imprints of the various modes will be discussed 
in Sec. VI. 


B. The Stress-Energy Pseudo Tensor for Gravitational Waves 


For all Lagrangian-based theories a very general method, with roots 
going back to Noether, 50 exists for constructing conserved quantities (see 
Ref. 5 and the references quoted therein for a more complete discussion). 
Invariance of the gravitational Lagrangian under coordinate transformations 
leads to the following identities: 


V IT =0 , 


ft' v - u v_ ) 

V|i U pA tG , v 

A . 

where ^ is the gravitational Lagrangian density, is the variational 

derivative of with respect to field y A occurring in ^ G , 


(76) 


•' v = 
' H 


v 


(77a) 


arid U V is defined by the functional changes of the y A , 8y A , under infini- 
(jA 

tesimal coordinate transformations, i.e., 

= x 11 + f , 


sy A = V V ,v - y A,/ 


(77b) 

(77c) 
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We have assumed ^ contains no higher than first derivatives of the w 
generalization to higher derivatives is straightforward. Equations (76) 
are of the form of conservation laws and our object is to identify in a 
physically meaningful way the gravitational portion of the conserved quantity. 
To facilitate the computation, we assume s^, has been rewritten in terms of 
and g^ [which can be done in principle by solving for h^g^, q^)]- 
Using the tensor transformation law for g^ and q^, one easily shows 


V V = - %(Qpp) V f ° r y A = *0*3 » 


(78) 


where parentheses denote symmetrization of indices. Using Eq. (78) we 
find the relation 

VV = ‘ 2 V(c£e) V K/ 5 W ' 2 V(c&s) V K /88 of> ) > (79) 


where 


8?q(ti, |) = 


If we now use the field equations 


^/ 88 og = - af^/ug^ 


(80) 


and Eq. (Uc), Eq. (79) becomes 


v „,A 

U . ^ ’ = 

|jA G 


2 V(cfe) 


V (^/ S n oe ) + (- g ) 1/2 t„ v * A , V + <- 8 > V V • < 81 > 
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We point out that although Eqs. (76) are "strong conservation laws"' 
(identities), one must use Eqs. (80) to get out a physically useful result. 
Substitution of Eq. (8l) into Eq. (76) yields 


(t V - (- g) 1//2 T V ) = 0 , 

' \i \i /, v 7 


(82a) 


3b 



where 


t' 1 




(82b) 


The conserved energy momentum vector is then 


l-i 


; (v° - <- ^\°) * 


(83) 


Since P in Eq. (83) contains a contribution from the matter stress energy 

U n 

tensor, we know we are on the right track. Problems arise when we notice 
that the quantity defined in Eq. (82b) is in general not positive definite, 
as a result of contributions from /\^ V . However, it can be shown from the 
generalized Bianchi identities of this theory (see Appendix B) that y\^ 
obeys the equation 


'V |V 


= 0 . 


(8M 
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Actually, Eddington was the first to point out that conservation laws of 


the form of Eq. (8^) follow from theories with absolute objects. If we 
now choose to work in the coordinate system in which is the globally 
constant Minkowski metric, Eq. (8h) becomes 


,■ 


- 0 , 


(85) 


and we see that y\^ V is conserved by itself, independently of energy gain or 
loss from matter (T V ). Since our usual idea of total energy conservation 
involves interactions, it is perhaps more useful to omit the separately con- 
served from consideration and to define, in this frame, the gravitational 

stress-energy tensor as 


t V = t ' v 


( 86 ) 


Thus y\^ V represents the energy density of a quantity associated with the 
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absolute field q ; at present we must regard it as a purely mathematical 

quantity whose noninteraction with matter mirrors the absolute nature of 

n . (As an aside, there always exists a t V which is a real tensor and 
uv l- 1 

. . 5s 

not a pseudo- tensor in prior-geometric theories of gravity, ) 

We point out that in the linearized approximation Eq. (85) is always 
the expression of Eq. (8^) in all frames related to the global Minkowski 
frame by infinitesimal coordinate (gauge) transformations. We proceed by 
explicitly calculating t V for the linearized theory. From Eq. (77a) 


* - + ^ = 


Inverting the linearized relation between g^ and h^ [cf. Eq. ( 70) ] and 
taking the required partial derivatives, we find 


*Vs, l/ ds 0 p,v - + 2K(1 - 8Krl 


Using Eqs. (87), (88), and Eq. (2) for s^, we finally obtain 

t/ - ( 16w) — 1 ♦ fh>\ a ) - (89) 

Since h transforms as a tensor, the above expression is gauge invariant. 

|1V 

Equation (89) expresses a naturally defined stress-energy complex for the 
gravitational field. 

Consider the energy density in a plane gravitational wave 

h TO = A r0 e lk< ^ ; k a k“ = 0 . (90) 

Then the first two terms in Eq. (89) do not contribute to t V and one obtains 


t V oc k k V , 
U M. 


(91a) 
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with 


t 0 °= (8,r 1 [a(h op)0 ) 2 + £(h >0 ) 2 ] 


(91b) 


With the suggested values for a and f [cf. Eq. (30)], Eq 0 (91b) indicates a 
positive definite energy density. It is encouraging to note that for pure 
spin 2 waves (only ^ present), Eq. (91b) becomes, for a = \/b [cf. Eq. (30)], 


0 

; o . 0 

spin 2 


*N‘‘(Vo )2 ’ ( 16 * ) ' 1 ( V ° )2 


(92) 


which is identical to the corresponding expression in GRT. 


VI. THE GRAVITATIONAL CONSTANT AND FURTHER EXPERIMENTAL TESTS 
A. A Time- Dependent Gravitational Constant 

As discussed in Sec. II, a number of existing solar system experiments 
place upper limits on the cosmological boundary values of h^ v [cf. Eqs. 

(31)- (32)]. These constraints can always be satisfied in a given epoch. A 

more relevant point is the time dependence of the u , which is directly 
related to the time dependence of the gravitational constant G. With the 
choice of adjustable constants given in Eq. (30), and using the explicit 
functional forms for D Q , D, one finds from Eq. (27) and Appendix A that 

1 - Jq (19 Wi + 7w Q ) + 0(co 2 ) = G . (93a) 

Thus 

(-i) dG/dt « - Yg ( 19u> 1 /dt + 7du 0 /dt) . (93b) 

Shapiro et al. 33 have placed limits on the time dependence of the gravita- 
tional constant by comparing the periods of planets with the ticking rates 
of atomic clocks. Ihey find 
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|(-i) (dG/dt) | < 4 x 10" 10 /year . (94) 

This constitutes an experimental constraint on the magnitude of the time 

derivatives of go occurring in Eq. (93b). Preliminary results from our 

p. v 

cosmological solution 19 indicate that the time dependences of and 
satisfy Eq. (94), but an improved Shapiro experiment might still prove to 
be a crucial experimental test of our theory. 

B. Gravitational-Wave Experiments 

The analysis of the preceding section reveals two crucial new experi- 
mental tests of our theory involving gravitational waves — two tests which 
have blossomed from our current program (see introductory remarks in Sec. 

I) — two tests which emphasize gravitational wave detection as a powerful 

3 b 

new tool for probing metric theories of gravity. ’ The two tests are (i) 
time delay between simultaneously emitted gravitational and electromagnetic 
waves and (ii) polarizations of gravitational waves. 

Since gravitational waves travel along geodesics of the "fast metric" 

ti and electromagnetic waves travel along geodesics of the "slow metric" 

0f3 

e there should be a time delay in reception of the two waves — emitted, 
for example, in simultaneous bursts by a supernova explosion. For waves 
emitted at the center of the galaxy, an order of magnitude estimate indicates 

Time Delay ~ (m/ r ) galaxy * travel time) 

~ (5 X 10 -7 ) * (3 X 10^ light years) s= 5 days . (95) 

Much longer delay times would hold for the Virgo Cluster. 

Polarization information is also a crucial experimental test. Equations 
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(75) indicate a purely longitudinal mode (tg), mixed longitudinal-transverse 
quadrupole type modes (^j ty^)* a P ure ^y transverse "breathing mode ($ 22 ^* 
and the familiar transverse quadrupole modes of GRT If an observer 

knows the direction of the wave, he can use Eqs. (75) to unambiguously catalogue 
the modes. If he does not know the direction of the source, he can still draw 
some conclusions. For example, if displacements do occur in more than one 
plane, then either the longitudinal- transverse modes (t 3 > tj) are present, 
or the purely longitudinal mode (\|Tg) is mixed in with one of the purely trans- 
verse modes (\|/^, *22^" 

It is important to note that until the problem of the generation of the 
various types of waves by particular sources is solved, our theory can only 
be verified by the presence of - but not ruled out by the absence of - the 
various possible modes indicated in Eqs. (75). This is unfortunate. But 
new doorways have been opened in the area of experimental tests and it is 
clear that gravitational tests outside of the PPN formalism must be contem- 
plated in the future. 
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APPENDIX A 


CONSTANTS APPEARING IN PN LIMIT (Sec. Il) 


“o = w oo 

W 1 s U 11 

CO = 3U)^ - Wq 


Eq. (12a): 


Eq. (12b): 


Eq. (12c): 
Eq. (13): 


D. = 


1 - 2Ko> + K 2 ^ 2 + 2Kww 0 + -| <o Q 2 - <o Q 


E- = 


1 - 2Kw - 2 w o 


F H - 2K + 2K to + 2K to r 


D 

E 

F 


2 2 3 2 

1 - 2Klo + to^ + K to - 2 Kioto ^ to^ 

1 * 2Kto + | ^ 


s - 2 k( 1 + to^ + 2K‘ 


00 


H = 1 - 2Kto - <o Q + to x 


^ D y 2 d- 3 / 2 


I, = 


I„ = 


1 /E M 

2 \D + D 0 / 

1 /3F fo E\ 

2 \ D " D Q + D/ 



L = - (a + Uf) -1 [f(l - 2Kto) + 2Ka(l - Kto)] 
M = - (a + 4f) _1 (2Ka + § f) 

N = 2k(f + Ka) (a + 4f) -1 
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Eq. (I6d) : 


S Q = 1^1 - 2Kco + L - § « 0 - M Wq ) - | - M 

S = I (1 - 2Ku + L - § u 0 - Mco q ) + N - 2K 

B q s I 3 (l - 2Kw + L - | cj q - Mu 0 ) - L - Mo^ 

. . 3 

Eq. (I6e): R Q = 1 - 2Kio + LJ 1 

Rj- = I 1 (Mw q - L) + M 

R 2 = Ig(Mu) 0 - L) - N 

= I_(Mu> 0 - L) + L + Mo^ 

Eq. (20a): Kj. = E^ - F Q (3C 1 - C Q ) 

K 2 = - [^(SCj - C 0 ) 2 + 2KC 0 (3C 1 - C 0 ) V f C 0 2 ] 

h s T[s o c o + S 1 (3C X - V 1 (E o + V - 3 tF o' rc o 
\ S t ^ E 0 B 0 ■ V R o + 3B 1 " B 0>1 

Eq. (20b): Kg s ECj^ + F(3C 1 - C Q ) 



APPENDIX B 


RELATIONS FOLLOWING FROM GENERALIZED BIANCHI IDENTITIES 


Assume that L„ has been rewritten as a function of and g^. Since 


P-v 


L is a scalar, its variation under infinitesimal coordinate transformations 
G 


must vanish, i.e., 


/ Ssfp \ Ij. 




l op " T ' 

Under the coordinate transformation 


-a a t a 

x = x + | , 


the functional changes in the tensors rj and g are 

= 


‘ ,p " ,a ’ 


2 T (a|P), 


S8qP = " 21 


Now define 


. - i (- I) l/2 T* • 


(Bl) 


(B2) 


. = -P 

a n ap § ’ 

(B3) 

; (a;P) 

(B4) 


(B5) 


<~G' u, 'ce 

and use the field equations to write 
Using Gqs. (B2)-(B6), Eq. (Bl) can be written in the form 

0 = 4- q) 1/2 P°\,V - 1 (- 8 > 1/2 < tC Vp + i ( - 8)1/2 5 

- (- u) 1/2 t 0 

Now if we remember that 


(B6) 


d\ 


(B7) 



y 


(B8) 


(^aV = (_ ,,) " 1/2 


(- ^ r^l 


a 


,3 


and also the corresponding equation for the covariant derivative with respect 

to g , the first two terms in (B7) vanish with proper boundary conditions 
6 Qp ? 

on Now use the matter equations, Eqs. (5), and the arbitrariness of 


l (and hence x ) to get from Eq. (B7) 

op 


IP 


= 0 


(B9) 


Equation (B9) is not an identity; we had to use both the matter and gravi- 
tational field equations to obtain it. [We would have obtained an identity 
in the place of Eq. (B9) had we not enforced the dynamical equations.] Since 
T) is covariantly constant with respect to "slash," Eqs. (B5) and (B9) imply 
the desired relation 


•Ajj. |v - t 2 t V(c£ p)( 5a TG // ' ST1 cp^ |v 0 


(BIO) 
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TABLE I 


Comparison of Construction of Stellar Models 


GRT Two-Metric 

Theory 


Number of coupled differential equations 
which must be integrated to find star's surface 


Type of differential equations used in deter- 
mining metric functions 


First-order 

linear 


Second order 
nonlinear 


Number of quantities whose central values must 
be chosen to satisfy boundary conditions 


Analytic Solutions 


Yes Probably 

not 


Uniqueness of solution for given central pres- 
sure and equation of state 

Number of parameters in exterior metric 
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FIGURE CAPTIONS 


F xg. 1 • 


Fig. 2. 


Fig. 3. 


Fig. U. 


Fig. 5. 


Radial Geodesics of Particles and Photons. The number along the 
curves indicate proper time values for massive particles released 
from rest at p = 10, 5. One of the curves is a photon geodesic and 
all curves have t -> oo and affine parameter -»■ w as p ■+ 1.5. 

Proper Surface Area of Sphere p = Const. The upper abcissa gives 
the proper time of an observer released from p = 10 as a local 
coordinate marker. 

Embedding Diagram for Equatorial Geometry. Solid line indicates 
Euclidean embedding (refers to z ordinate) and dashed line indi- 
cates pseudo-Euclidean embedding (refers to iz ordirate). Numbers 
along curve indicate values of p. 

Effective Potential for Massive Objects. Dots indicate circular 
orbits. 

Effective Potential for Photons. 
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